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Abstract 

We investigate the superfluid phase transition in a gas of Fermi atoms loaded on a three- 
dimensional optical lattice. When the lattice potential is strong, this system can be well described 
by an attractive Hubbard model. In this model, we calculate the supernuid phase transition tem- 
perature T c , including both supernuid and (spin and charge) density fluctuations within the self- 
consistent i-matrix theory and fluctuation exchange approximation, respectively. Since we treat 
these fluctuations in a consistent manner, our theory satisfies the required particle-hole symme- 
try over the entire BCS-BEC crossover region. We show that charge density fluctuations compete 
against supernuid fluctuations near the half-filling, leading to the suppression of T c . As a result, the 
maximum T c is obtained away the half-filling. Since the strong density fluctuations originate from 
the nesting property of the Fermi surface at the half filling (which is absent in a uniform gas with 
no lattice potential), our results would be useful in considering lattice effects on strong-coupling 
superfluidity. 

PACS numbers: 71.10.Ca, 03.75. Ss, 37.10.Jk 
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I. INTRODUCTION 



Optical lattice is an artificial lattice produced by standing wave of laser light [1]. In a cold 
atom gas loaded on an optical lattice, atoms feel a periodic potential due to the Stark effect. 
When the lattice potential is strong, this system can be well described by the Hubbard 
model, where atoms are hopping between nearest-neighbor sites, interacting with each other 
when they meet at the same lattice site. Since the Hubbard model is a fundamental model 
in condensed matter physics, it is expected that various topics discussed in this field may be 
solved by using the optical lattice system. Indeed, the superfluid-Mott insulator transition 
has been observed in a 87 Rb lattice Bose gasfl [3]. More recently, the superfluid state has 
been also realized in a 6 Li lattice Fermi gasjjil. 

Besides the optical lattice, a tunable interaction associated with a Feshbach resonance is 



also an advantage of cold atom gases 6|, Q| • Using this unique property, several experimental 
groups . fioi. Tll| have succeeded in realizing Fermi superfluids and the BCS-BEC (Bose- 

Einstein condensation) crossover 0, Q, Q, Q, 0, H 0, Q 

in the absence of optical 

lattice. The BCS-BEC crossover is a very interesting phenomenon, because one can study 
the weak-coupling BCS-state and the BEC of tightly bound molecules in a unified manner by 
varying the strength of a pairing interaction. Since this tunable pairing interaction also works 
in an optical lattice, it is interesting to examine how the BCS-BEC crossover phenomenon 
is observed in the optical lattice system. We briefly note that, in the Hubbard model, 
interaction effects are parametrized by the scaled interaction U/t (where U and t represent 
an on-site pairing interaction and nearest-neighbor hopping, respectively). Thus, in addition 
to the direct tuning of the pairing interaction by a Feshbach resonance, continuous change 
from the weak- to strong-coupling regime can be also realized by adjusting the hopping 
parameter t by tuning the intensity of laser light producing the optical lattice. 

The BCS-BEC crossover in the attractive Hubbard model has been discussed in super- 
conductivity literature [3] , in connection to strongly-correlated electron systems. Nozieres 
and Schmitt-Rinkjisl pointed out the mass enhancement of tightly bound molecules in the 
strong-coupling BEC regime, because of virtual dissociation of a bound molecule during hop- 
ping between lattice sites. This mass enhancement is expected to decrease the superfluid 
phase transition temperature T c in the BEC regime, which has been theoretically confirmed 
by the self-consistent t-matrix theory [221] . dynamical mean-field theory [23], and quantum 
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wave (CDW) fluctuations 



Monte-Carlo simulation 24 L Competition between pairing fluctuations and charge density 



251 ] near the half-filling (which comes from the nestin 



of the square-shape 



dimension, Refs. 



27 



ermi surface at the half-filling) has been also studied 



g property 



26t ] . In two- 



28] pointed out that this competition leads to vanishing T c at the 
half-filling. For more details, we refer to Ref. 

In this paper, we investigate the superfluid phase transition in a gas of Fermi atoms 
loaded on a three-dimensional cubic optical lattice. Treating this system as the Hubbard 
model, we calculate T c in the BCS-BEC crossover region, including pairing fluctuations, as 
well as CDW and spin density wave (SDW) fluctuations in a consistent manner within the 
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29( | and fluctuation exchange approx- 



self-consistent i-matrix approximation (SCTA) 
imation (FLEX)|28j, respectively. Our theory satisfies the required particle- hole symmetry 
over the entire BCS-BEC crossover region. While a finite T c is obtained at the half-filling 
in contrast to the two-dimensional case 27J] , the superfluid phase transition is shown to be 



strongly influenced by CDW fluctuations near the half-filling. The resulting T c takes the 
maximum value, not at the half-filling, but around the quarter filling. Since strong CDW 
fluctuations are characteristic of the lattice system we consider in this paper, the observation 
of the filling dependence of T c would be an interesting problem. 

This paper is organized as follows. In Sec. II, we explain our formulation. The self- 
consistent t-matrix approximation (SCTA) for pairing fluctuations and the fluctuation ex- 
change approximation (FLEX) for CDW and SDW fluctuations are explained. In Sec. Ill, 
we present our numerical results for the superfluid phase transition temperature, only taking 
into account pairing fluctuations. We examine effects of CDW and SDW fluctuations on the 
superfluid phase transition in Sec. IV. Throughout this paper, we set % = ks = 1- 

II. FORMULATION 



We consider a two-component Fermi gas in a three-dimensional cubic optical lattice. In 
superfluid Fermi gases, all the current experiments are using a broad Feshbach resonance 0, 



S0 



In this case, details of the Feshbach resonance is known to be not important as 



far as we consider the interesting BCS-BEC crossover region, so that we can safely consider 
this system using the ordinary BCS model. In addition, as mentioned in the introduction, 
a Fermi gas in an optical lattice can be well described by the Hubbard model when the 
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FIG. 1: Self-energy corrections describing (a) pairing fluctuations, (b) CDW fluctuations and SDW 
fluctuations of the z-component, and (c) SDW fluctuations of the x- and y-component. The solid 
line and dashed line describe the single-particle Green's function G and the attractive interaction 
—U, respectively. 
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FIG. 2: Vertex functions of (a) particle-particle channel r pp and (b), (c) particle-hole channel 



T\. In our calculation, the lowest order term in terms of U is included in T 



pp • 



lattice potential is strong. Under these conditions, we consider the attractive Hubbard 
model described by the Hamiltonian 



h = -t ( c l* c j,° + h - c ~ U J2 n ^ n n - ^ Yl 

(i,j),<r i i,er 



(1) 



Here, c JjCr is the annihilation operator of a Fermi atom at the j-th lattice site, where the 
pseudo-spin o =|, \ describes two atomic hyperfine states, t is the hopping matrix element 
between nearest-neighbor sites, and the summation in the first term is taken over 

nearest-neighbor pairs. rii a = c^c^ is the number operator. The on-site pairing interaction 
—U (< 0) is implicitly assumed to be tunable by using a Feshbach resonance. \i is the Fermi 
chemical potential. In Eq. (pQ), we have neglected effects of a harmonic trap, for simplicity. 



We calculate the superffuid phase transition temperature T c , extending SCTA developed 
by Haussmannjl6. fir! ] so as to include CDW and SDW fluctuations within FLEX. Effects of 
these fluctuations are described by the self-energy T,(k,iu m ) in the single-particle Green's 
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function G(k,iu m ), 

G(k, iu m ) = 1 (2) 

where £k = —2t\ cos /c x + cos hy -\- cos fc z ] — fi is the kinetic energy of a Fermi atom measured 
from the chemical potential \x (where the lattice constant is taken to be unity). uj m is the 
fermion Matsubara frequency. 

Figure Ufa) shows the self-energy correction coming from pairing fluctuations (= S pp ). 
Summing up this type of diagrams within SCTA, we obtain 

S pp (fc, iuj m ) = 4 H r PP (qr, iu n )G(q - k, iv n - iu m ), (3) 

P q,iv„ 

where v n is the boson Matsubara frequency. Hereafter, the factor 1/N of lattice sites number 
in front of momentum summation is abbreviated for simple presentation. r pp is the particle- 
particle scattering vertex diagrammatically described by Figj2]^a). The result is 

r » (9> " ) = -i-ra! (g >„r (4) 

where, 

n pp (q, iu n ) = - G ( k i iu m )G{q - k, iv n - iuj m ) (5) 

P k,«w m 

is a correlation function describing fluctuations in the Cooper channel. 

In Fig{TJ panels (b) and (c) describe fluctuations in the particle-hole channel within FLEX. 
Panel (b) involves both CDW fluctuations and longitudinal SDW fluctuations. Panel (c) 
involves transverse SDW fluctuations. Summing up these diagrams, we obtain the self-energy 
corrections associated with CDW fluctuations (= £ ph ) and SDW fluctuations (= S ph ) as 

£ p h(fc, iujm) = - — Yl n ph(<?> ^n)r p h(q, iv n )G{k - q, iu m - iv n ) 

P q,iv n 

= \ V ph(<l> i"n)G{k - q, iu m - iv n ), (6) 

where V^ s is introduced for saving the computational time by performing the fast Fourier 
transformation (FFT), which is explained in Sec. III. Here, the vertex functions T^ h and r ph 
are obtained from the sum of the diagrams shown in Figsf2]^a) and[2^b), and their expressions 
are given by 
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F s (n iu 1 - 3 ^Ph^^n) (R s 

T ^ ZUn) -2i + uu ph ( q ,i, n y (8) 

In Eqs. and (JBJ), the correlation function n p h(fe, iu n ) describes fluctuations in the particle- 
hole channel, having the form 

Il ph (qr, iv n ) = iuJ m)G(k - q, iu m - iv n ). (9) 

' k,iui m 



The superfluid phase transition temperature T c is determined from the Thouless criterion la . 



171 ]. stating that the superfluid phase transition occurs when the particle-particle scattering 



vertex r pp (q,iz/ n ) has a pole at q = v n = 0. Using this, we obtain the equation for T c as 

l = UIL pp (q = 0,iv n = 0). (10) 

We note that Eq. (JTUI) is affected by CDW and SDW fluctuations through the self-energy 
£ = S pp + Sp h + S ph in the Green's function. 

In the weak-coupling BCS regime, we may set \i = e-p (where e-p is the Fermi energy) in 
Eq. ( JTUI) . However, the chemical potential is known to deviate from e-p, as one approaches 
the strong-coupling BEC regime 12, ll3|]. This strong-coupling effect is taken into account 
by considering the equation for the filling number n (which gives the number of atoms per 
lattice site), given by 

n =l E e^ s G(k,icu m ), (11) 

r k,iLU m 

where 5 is an infinitesimal positive number. We solve the coupled equations (JTUI) and (JTTj) 
to determine T c and \x self-consistently for a given U and n. 

At the half-filling n = 1, the superfluid state and CDW are degenerate in the sense that 



they have the same phase transition temperature |2CM. |26|]. The CDW phase transition is 
characterized by the divergence of the charge susceptibility Xcdw(Q) with the momentum 
Q = (7r, 7T, 7t). In the random phase approximation, Xcdw is given by 

(n ( \\ 1 n ph (Q,0) 
Xcdw(Q = (-, tt, 7T)) = - 2 1 _ uuM0 y (12) 

where the correlation function n p h(Q, iv n ) is given by Eq. (J9j). Comparing Eq. ( JT2l) with 
Eq. ((7j), we find that the CDW vertex function Tp h (Q, zz/ n = 0) also diverges at T c . This 
CDW instability at the superfluid phase transition temperature T c is absent when 1 due 
to the absence of the perfect nesting of the Fermi surface. However, since the denominator 
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1 — UHp^QjiVn = 0) in Eq. (fT2l) is still small at T c near the half-filling, strong CDW 
fluctuations are expected when n ~ 1. 

In contrast, the SDW vertex function T ph in Eq. does not diverge at T c even when 
n — 1. Namely, spin fluctuations are weak in the attractive Hubbard model. 

III. EFFECTS OF PAIRING FLUCTUATIONS ON T c AND fi IN THE BCS-BEC 
CROSSOVER REGION 

In the following two sections, we show our numerical results obtained by solving the 
coupled equations (fTUj) and (TTTT) . In this section, we first consider T c in the BCS-BEC 
crossover, including pairing fluctuations only. Although one cannot actually ignore strong 
CDW fluctuations near the half-filling, examining this simple case is still useful in considering 
importance of CDW and SDW fluctuations. We separately discuss effects of CDW and SDW 
fluctuations in Sec. IV. 

Before showing our results, we summarize the outline of computation. In solving the 
coupled equations (fTOl) and ffTTj) , we use the fact that the self-energy S = £ pp + Sp h + S ph in 
Eqs.Q and the correlation functions n pp and II ph in Eqs. ([5]) and ([9]), and the number 
equation in Eq. ( TTTT) have simple expressions in real space, as 



r pp (r, r)G(-r, -r) + V d h (r, r) + VUr, r) G(r, r), (13) 



S(r,r) = S pp (r,r) + Sj h (r,r) + S ph (r,r 



U pp (r,r)=G(r,r)G(r,r), (14) 

n ph (r,r) = -G(r,r)G(-r,-r), (15) 

n = 2G(r = 0, r = -5). (16) 

Here, r is the spatial position of a lattice site and r is the imaginary time. The Fourier 
transformation is defined by 

G ^ r ^ = \ £ G(k,iu m )e^ r -^ T \ 

p (17) 

—i(k-1 — u) m -r) 



Ef 
/ drG(k,T)e 
r Jo 
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FIG. 3: Superfluid phase transition temperature T c as a function of pairing interaction U and 
filling number n. In this figure, and in FigUl we only include pairing fluctuations described by 
E pp . Since T c at the filling 2 — n is the same as T c at n due to the particle- hole symmetry, we only 
show the result less than half-filling (n < 1). 

I — I 

To use Eqs. ( IT3l) -( fT6l) . we employ the FFT method [301] . We discretize the momentum region 
< k x , k y , k z < tt into 16 x 16 x 16 cells. For the frequency summations, we introduce a 
finite cutoff frequency cj max = 7rT(2n max + 1) for fermions and z/ max = 2iiTn Tnax for boson, 
with n max = 512. The values of these cutoffs are chosen so as to be much larger than the 
band width 2zt (where z = 6 is the coordination number of the simple cubic lattice), as well 
as the magnitude of the pairin g in teraction U . To avoid effects of these cutoff frequencies, 
we use the method discussed in [28[ . We explain the outline of this method in the Appendix. 

Figure [3] shows the calculated T c in the BCS-BEC crossover. In this calculation, we only 
include the self-energy S pp associated with pairing fluctuations. Since the Hubbard model 
has the particle-hole symmetry, the filling dependence of T c is symmetric with respect to 
n — 1 (although we do not show it explicitly). Namely, for a given U, Figj3] shows that the 
maximum T c is obtained at the half-filling n = 1, while T c vanishes at n = 2 as in the case 
of n = 0. 

In the weak-coupling BCS regime (U/Qt < 1), T c is an increasing function of U. This 
behavior agrees with the well-known mean- field BCS result, 

Tcocte" 1 ^ 7 , (18) 

where N(0) is the density of states at the Fermi level. In this regime, since T c is comparable 
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FIG. 4: Chemical potential \i as a function of U at T c . u = —U/2 is the exact solution at the 
half-filling n = 1. At U = 0, apart from weak temperature effect, fi is almost equal to the Fermi 
energy ep for a given filling number n. 

to the binding energy E'bind of a Cooper pair at T = 0, the increase of T c reflects the 
enhancement of -Ebind as one approaches the strong-coupling regime. 

For a given filling number n, the maximum T c (= T c max ) is obtained around the inter- 
mediate coupling region U/6t ~ 1. At the half-filling (n = 1), we obtain T c max = 0.042 
at U/6t = 0.9. T c then decreases as one further increases the magnitude of U. The de- 
crease of T c in the BEC regime is characteristic of the BCS-BEC crossover in the Hubbard 



model 20|]- In a uniform Fermi gas with no lattice potential, T c approaches the constant 
value T c = 0.218e F in the BEC limit [l5, 3, Q- 

We note that the decrease of T c in the BEC regime does not mean the small binding 
energy of a bound molecule in this regime. As shown in FigJU while the chemical potential 
H is almost equal to the Fermi energy e-p at U = (apart from weak temperature effect), 
it becomes smaller than the bottom of the band (= —Qt) in the strong coupling regime. 
This means the existence of a finite energy gap E g = |//| — 6t in the Fermi single-particle 
excitations in the BEC regime. Since E g is directly related to the dissociation energy of a 
bound molecule, we find that the binding energy continues to increase even in the BEC 
regime (although T c decreases as shown in Figj3]). 

To understand physics behind the decreasing T c in the BEC regime, it is helpful to derive 
an effective model valid for this regime. To do this, we note that in the BEC regime, 
tightly bound molecules have been already formed above T c . In this pointed out in 
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Ref . [13j , molecular motion is accompanied by virtual dissociation, because each atoms in a 

molecule has to move one by one in the Hubbard model. In addition, this virtual dissociation 

I 

also leads to a repulsive interaction between molecules 13]. Including these effects within the 



second order perturbation in terms of the hopping t, we obtain the effective Hamiltonian 



2l|, 



tfeff = ~ £ l4bj + h.c] + -^£ nfn] -MbE nf, (19) 

where b\ = cjfcjj describes a molecule at the z-th lattice site, and nf = (n^ + nii)/2 gives the 
number of molecules under the assumption that all the atoms form small on-site bound pairs. 
/ib = 2{i + U+2zt 2 /U is the molecular chemical potential to control the molecular density. In 
Eq. f[T§j) . double occupancy of molecules is forbidden due to the Pauli's exclusion principle 
of Fermi atoms in them. Noting this and commutation relations, = 2(nf — 1/2), 

[nf - 1/2, &J] = b\, and [nf - 1/2, & 4 ] = -b h we can map Eq. (JTSD onto the 5 = 1/2 
Heisenberg model, by replacing (b\,bi,nf — 1/2) with ((— l) 4 ^, (—iyS^,Sf), 

H eS = J^S i -S j -h B Y,S*. (20) 

Here, J = 4t 2 /U is an exchange interaction and Hb = /^b — 2zt 2 /U works as an external 
magnetic field. Since Sf = —1/2 and Sf = +1/2 state in Eq. (Tl6|) . respectively, correspond 
to a vacant and occupied site in the effective model in Eq. (Tl9l) . the half-filling case is 
described by setting h-& = in Eq. (|20|) . Equation (1201) clearly shows that the 'Neel 
temperature Tn' (which corresponds to T c in the original Hubbard model) is lower for larger 
U, consistent with the decreasing T c in the strong-coupling BEC regime shown in FigfSJ 

The Neel temperature of the Heisenberg model has been studied by various methods. It 
has been shown that the mean-field result T N = Qt 2 /U is suppressed by spin fluctuations to 



be Tn = 3.59t 2 /U (high temperature expansion) [311] and Tn = 3.78t 2 /U (quantum Monte 



Carlo method) 32|. In the present calculation, we obtain T c ~ 1.5t 2 /U < (3.59 ~ 3.78)t 2 /U 
in the BEC regime, showing that the present calculation based on SCTA underestimates T c 
in the BEC regime. 

One reason for this discrepancy is the approximate treatment of the interaction between 
molecules. To see this, we consider the BCS-BEC crossov er p roblem within the Gaussian 



fluctuation theory developed by Nozieres and Schmitt-Rink 



131 ] . (We refer this theory as the 



NSR theory in the following.) In the NSR theory, one also solves the coupled equations (EG) 
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FIG. 5: Calculated T c within the NSR theory. 'BCS limit' is the mean-filed BCS result with 
H = £p. 'BEC limit' shows T c of an ideal molecular Bose gas. 

and (TTTj) . where the single-particle Green's function G and the correlation function n pp are 
now replaced by 

G(p, iu m ) = G (k, iu m ) + G (k, iu m )T, pp (k, iu m )G (k, iu) m ), (21) 
U pp (q,iu n ) = - J2 G (k,iu; m )G (q-k,iiy n -iuj m ), (22) 

r k,iw m 

where Go(k, iu m )~ l = iojm—e^+n is the Green's function in a free Fermi gas. The self-energy 
part in Eq. f[2~Tj) is given by 

S pp (fc, iu m ) = - F pp(Q> iv n )G (q - k, iv n - iuj m ). (23) 

The particle-particle scattering vertex r pp is given by Eq. (j4j), where n pp is replaced by 
Eq.fl22D. 

Figure \5\ shows the calculated T c based on the NSR theory. Comparing this result with 
Fig 01 we find that the NSR theory gives higher T c . Since the NSR theory is a low den- 



sity approximation [l3j (Note that only the free propagator Go is used and the self-energy 
correction is only taken into account to the first order in the number equation (jTU]) .). the 
difference between the two is more remarkable in higher filling cases. 

In the BEC limit of the NSR theory, Eq. (|T0l) gives fi = —U/2, and the number equation 
( TTTT) reduces to the condition for BEC in an ideal Bose gas, 

2 = ? ^Ffe) _ l' ^ 
11 



where E® = —2(2t 2 /U)[cosq x + cosq y + cosg 2 ] and p,B = —6(2t 2 /U). Indeed, Fig J5] shows 
that T c obtained from Eq. (1241) well describes the NSR result in the strong-coupling regime. 
Noting that the kinetic energy can be also obtained from the first term in Eq. (TT9|) 
when one regards hi as a boson operator, we find that the NSR theory ignores the repulsive 
interaction between molecules given by the second term in Eq. (1191) . 

In a uniform Fermi gas with no lattice potential, Haussmann pointed out that, in 
the BEC regime, SCTA includes the interaction between molecules within the Born 
approximation 16|. This molecular interaction can be written as U-q = Ana-B/M^, where 



Mb is a molecular mass. The s-wave molecular scattering length a-Q is related to the s-wave 
atomic scattering length a s as as = 2a s . Even in the presence of the lattice, the molecu- 
lar interaction is expected to be included within the same approximation level. Thus, we 
find that the lower T c within SCTA than the NSR result originates from the molecular 
interaction. Namely, the molecular interaction lowers T c in the BEC regime in the lattice 
system. 



Recent work 



33 



34 



35 



361 ] on a uniform Fermi gas has clarified that, when one carefully 



treats higher order molecular scattering processes and a finite value of molecular binding 
energy, a-Q reduces to = O.Qa s < 2a s . This clearly indicates overestimate of the magnitude 
of molecular interaction in SCTA. When we apply this discussion to the present lattice 
system, one reason for the underestimate of T c (= IM 2 /U < (3.59 ~ 3.78)t 2 /U) is expected 
to be the overestimate of the molecular interaction. When one could correct this point, T c 
would be higher to be close to the 'Neel temperature' of the Heisenberg model in Eq. ( |20l ). 
This improvement is an interesting problem; however, in this paper, leaving this as a future 
problem, we treat pairing fluctuations within SCTA and discuss effects of CDW and SDW 
fluctuations in the next section. 

Before ending this section, we briefly note that the NSR theory does not satisfy the 
particle-hole symmetry, when it is applied to the Hubbard model. As shown in FigJHJ the 
calculated T c based on the NSR theory is unphysical around n = 2. (Note that the system 
must be a band insulator at n = 2, leading to vanishing T c .). We also find that the required 
symmetric filling dependence of T c with respect to n — 1 is not obtained within the NSR 
theory. In addition, although the chemical potential must satisfy /i = —U/2 at n = 1 due 



to the particle-hole symmetry 



371 ]. the NSR result satisfies it only in the BEC limit U — > oo, 



as shown in FigJTl In contrast, in addition to the symmetric filling dependence of T c , SCTA 
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FIG. 6: T c as a function of the filling number n. The result within the SCTA and that within the 
NSR theory are compared. We set U/6t = 1. 
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FIG. 7: Chemical potential \x at T c calculated within the NSR theory. 



can also reproduce the exact result [i = —U/2 at n = 1 over the entire BCS-BEC crossover. 
(See FigHl) We emphasize that satisfying these required conditions is important in any 
consistent theory. 



IV. EFFECTS OF CDW AND SDW FLUCTUATIONS 

Figures [8] and [9] show calculated T c and /x(T c ) respectively in the BCS-BEC crossover, 
when the CDW (Sp h ) and SDW fluctuations (S ph ) are both taken into account. As shown 
in Fig. [9l the required condition fi = —U/2 at n — 1 is still satisfied when one includes 
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FIG. 8: Superfluid phase transition temperature T c as a function of pairing interaction U and 
filling number n. In this figure, and in FigJH we include CDW and SDW fluctuations described by 
in addition to pairing fluctuations. In comparison with FigHJ T c around n = 1 are suppressed 
in the strong coupling regime. 
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FIG. 9: Chemical potential // as a function of U at T c when we consider CDW and SDW fluctuations 
in addition to pairing fluctuations. Self-consistent treatment of CDW and SDW fluctuations does 
not break the particle- hole symmetry condition fj, = — U/2 at n = 1. 

CDW and SDW fluctuations within FLEX. 

Comparing FigJUwith Figj9j, one finds that effects of CDW and SDW fluctuations on the 
chemical potential \i are weak. In contrast, from the comparison of Figj3] and FigJHl T c is 
found to be suppressed near the half-filling when these fluctuations are taken into account. 

To see the suppression of T c more clearly, we show T c at n = 1 in FigJTUl At the half-filling, 
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FIG. 10: Effects of CDW and SDW fluctuations on T c at the half-filling n = 1. SCTA: Pair- 
ing fluctuations are only included. SCTA+CDW+SDW: Pairing, CDW, and SDW fluctuations 
are all included. SCTA+CDW: Pairing and CDW fluctuations are taken into account. These 
abbreviations are also used in Figs [IT] and [TZl 
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FIG. 11: Denominator of the static charge susceptibility x(Q) i n Eq. (|12p at the superfluid phase 
transition temperature T c . The CDW instability is determined when 1 — UU p h(q = Q,u = 0) = 
is realized. 

the CDW instability occurs simultaneously (which is confirmed by vanishing denominator of 
the charge susceptibility x(Q) as shown in Fig JTTl) . which leads to the remarkable suppression 
of T c . However, in contrast to the two-dimensional case, where T c vanishes at n = 1[27], we 
still obtain a finite T c even at the half-filling. 

Although SDW fluctuations are weak in the attractive Hubbard model, we still find 



SCTA(n=1.0) 

/"n SCTA+CDW+SDW(n=1.0) 

SCTA+CDW(n=1.0) 




SCTA+CDW+SDW(U/6t=0.7) 

SCTA+CDW+SDW(U/6t=2.0) 

SCTA+CDW(U/6t=0.7) x 
SCTA+CDW(U/6t=2.0) □ 
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FIG. 12: Filling dependence of the superfluid transition temperature T c . Upper and lower panels 
show the cases of the weak-coupling BCS regime (U/6t = 0.7) and strong-coupling BEC regime 
(U/6t = 2), respectively. 

their effects around U/6t ~ 1 in FigflQl From the comparison of the result referred to as 
'SCTA+CDW+SDW with 'SCTA+CDW in FigM we find that SDW fluctuations weaken 
the suppression of T c by CDW fluctuations. 

In this intermediate coupling regime, since the binding energy of a Cooper pair is not 
very strong, (pseudo)spin degrees of freedom still remains, which contribute to SDW fluc- 
tuations. As one approaches the strong-coupling regime, these spin degrees from free- 
dom disappear due to the formation of singlet pairs. Indeed, in FigflQl the two results, 
'SCTA+CDW+SDW and 'SCTA+CDW, give almost the same T c when U/Qt > 3. 

Although the degeneracy of the superfluid state and CDW is absent when n ^ 1, we can 
still expect strong influence of CDW near the half-filling due to strong enhancement of charge 




SCTA+CDW(U/6t=0.7) 




SCTA+CDW(U/6t=2.0) 
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susceptibility shown in Fig{HJ Indeed, as shown in FigJIU CDW fluctuations suppress T c 
near n = 1. (Compare 'SCTA' with 'SCTA+CDW+SDW in FigiH) The maximum 
T c is thus obtained, not at the half-filling, but away from the half-filling. Although SDW 
fluctuations enhance T c , the overall behavior is unchanged. (Compare 'SCTA+CDW+SDW 
with 'SCTA+CDW in FigJT21) Since spin degrees of freedom are almost absent in the strong- 
coupling regime, effects of SDW fluctuations are weaker in the lower panel (BEC regime) 
than the upper panel (BCS~crossover regime) in FigJT2l 

V. SUMMARY 

To summarize, we have discussed the superfluid phase transition in the BCS-BEC 
crossover regime of a two-component Fermi gas loaded on a three-dimensional optical lattice. 
Treating this system as the attractive Hubbard model, we calculated the superfluid phase 
transition temperature T c , including pairing fluctuations within the self-consistent t-matrix 
theory, as well as CDW and SDW fluctuations within the fluctuation exchange approxima- 
tion. We determined T c and the Fermi chemical potential /i self-consistently as functions of 
the pairing interaction U and filling number n in the BCS-BEC crossover region, by solving 
the equation for T c , together with the number equation. 

CDW fluctuations are strong near the half-filling due to the nesting property of the 
Fermi surface at n — 1. These strong fluctuations remarkably decrease the superfluid phase 
transition temperature T c around n = 1. As a result, the maximum T c is obtained, not at the 
half-filling, but away from the half-filling. We also showed that, although SDW fluctuations 
are weak in the attractive Hubbard model, they still affect T c to increase slightly in the 
intermediate coupling region U/6t ~ 1. 

We have also discussed the validity of SCTA. Our theory satisfies the required condi- 
tion associated with the particle-hole symmetry of the Hubbard model which is not derived 
from NSR theory. Therefore, it would be a good starting point to improve the BCS-BEC 
crossover theory in optical lattices. On the other hand, we showed that this approximation 
underestimates T c in the BEC regime. As a key to understand this, we pointed out the 
importance of a repulsive interaction between molecules. Since the Gaussian fluctuation 
theory (which completely ignores the molecular interaction at T c ) largely overestimates T c , 
we expect that SCTA overestimates effects of the molecular interaction. Indeed, the overes- 
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timate of the molecule interaction within SCTA has been pointed out[16|, |33j, |34j, |35|, [36|] in a 
uniform Fermi gas with no lattice potential. Inclusion of the correct value of the molecular 
interaction in the present theory is our future problem. 

So far, the superfluid Fermi gas in an optical lattice has been realized when the lattice 
potential is weak[4j. To really realize the Hubbard model in a cold gas of Fermi atoms, 
stronger optical lattice potentials must be used, which, however, inevitably leads to a small 
nearest-neighbor hopping t, as well as low T c . Thus, to realize Fermi superfluids in such a 
difficult situation, our results indicate that the filling number should be set to be away from 
n — 1 (to avoid the suppression of T c by CDW fluctuations) so that one can reach T c as easy 
as possible under a given experimental condition. Since the Hubbard model is a fundamental 
model in condensed matter physics, realization of Fermi superfluid in the Hubbard model 
produced by strong optical lattice potential would be great challenge in cold atom physics. 



Acknowledgments 

This work was supported by a Grant-in-Aid for Scientific Research of Priority Area 
"Physics of New Quantum Phase in Superclean Materials" from the Ministry of Educa- 
tion, Culture, Sports, Science and Technology of Japan (MEXT). The authors were also 
supported by MEXT (No. 19340099 (H. T. and K. M.), and 19540420 (Y. O.)). 



APPENDIX A: APPLICATION OF FFT TO FREQUENCY SUMMATIONS 

In this Appendix, we explain how to use FFT algorithm in transforming between the 
Matsubara frequency and the imaginary time. As an example, we consider the single- 
particle Green's function G here. However, the method explained in this appendix can be 
also used in calculating the self-energies as well as correlation functions. 

We introduce a Matsubara frequency cutoff u; max = 7rT(2n max + 1) and evaluate the 
Green's function by using the Fourier transformation. By introducing this cutoff frequency, 
the Fourier transformation from the Matsubara frequency into imaginary time can be rewrit- 
ten as (suppressing the variables k and r), 

^max 

G(r) = T G(iu m )e-^ r 

m=-n max -l 
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= 2T (Ke[G(iu m )] cosu m r + lm[G(iuj m )) sincj m r), (Al) 

where we have used the analytic property G(— iuj m ) = G*(iu m ), as well as the fact that G{r) 
is a real function. Replacing r by (3 — r in Eq. flAll) . we have 

^max 

G{(3 — t) = 2T (— Re[(7(za; m )] cosa; m T + lm[G(iuj m )] sin QJ m r). (A2) 

m=0 

Thus, one may only consider the region 0<r</?/2in executing the cosine (sine) Fourier 
transformation from the Matsubara frequency into imaginary time. In numerical calcula- 
tions, we divide the region < r < (3/2 into n max cells and use FFT method. 

When we calculate the inverse Fourier transformation, we meet the problem that the ex- 
pected high-frequency behavior G{k, iu m ) ~ l/iuj m is not obtained because of the introduced 
cutoff frequency u max . To avoid this problem, we rewrite the inverse Fourier transformation 
in the form, 

2n m 



G(ioo m ) = / dTG(r)e iuJmT = ]T / drG{r)e^ T . (A3) 



'%TG(r)e luJmT = 

3=1 JT '- 1 

Here, Tj = Arj, where At = /3/2n max . When we approximately write the Green's function 
in the region r = [t 3 -_i,Tj] as G(r) ~ G(r 3 _x) + (r — 7j_i) x (G(tj) — G(r ? _ 1 ))/Ar, we can 
execute the integrals in Eq. (|A3|) . The result is 

G(iu m ) = —[-G(r 2n _ =(3-5)- G(r = +5)} + - J—J-^fa) - G(r )} 

+ £ {2G(r,) - G(r i+1 ) - G(r i _ 1 )}e^ - {G(r 2rtma J - G(r wl )} , 

(A4) 

where 5 is a infinitesimal positive number. In Eq. (1A4I) . because G(/3 — 5) = — G(— 5), the 
first term remains finite due to the discontinuity of Green's function at r = 0, giving the 
expected high frequency behavior (~ l/iu m ). We apply FFT to calculate the second term 



in (1A4I) . We note that Eq. (1A4I) has been also derived in Ref. 



281 ] by integration by parts. 
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